Quantised charge pumping through multiple quantum dots 
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We study electron pumping through a system of barriers, whose heights are deformed adiabatically. 
We derive a simple formula for the pumped charge Q in terms of the total reflection and transmission 
amplitudes and phases. The pumped charge increases with the number of barriers (nt,) and shows an 
interesting step-like behaviour, with the steps appearing at integer values of Q. The pumped charge 
also tends towards quantisation with the increase of the amplitude of the time-varying potential. 
The value of the quantised pumped charge is shown to be correlated to the discontinuity of the 
reflection phase. 

PACS number: 73.23.Hk, 72.10Bg, 73.40Ei 



A parametric electron pump is a device which gener- 
ates a dc current at zero bias by cyclic deformations of the 
system parameters [1,2]. The parameters of the Hamil- 
tonian are slowly varied as a periodic function of time 
such that the Hamiltonian returns to itself at the end of 
one cycle, but charge has been pumped through the sys- 
tem. In the last few years, electron pumps consisting of 
small semi-conductor dots have received a great deal of 
experimental [3,4] and theoretical attention [5-9]. 

Of more recent interest are the quantum pumps in open 
dot systems, where quantum interference of the electronic 
wave-function rather than Coulomb blockade (CB) is ex- 
pected to play the major role. Such a pump has been 
fabricated by Switkes et al [10]. A scattering approach 
to such a parametric electron pump was pioneered by 
Brouwer [6] where the pumped current was related to 
parametric derivatives of the scattering matrix. Using 
this, several theoretical papers [11-13] have investigated 
the connection between resonant transmission and the 
pumped charge. In general, the pumped charge is not 
quantised; however, it has been shown that when the 
pumping contour encloses almost all of a resonance line, 
the charge pumped is almost quantised [12]. It has also 
been shown that inclusion of inter-electron interactions 
using the Luttinger liquid formalism [14] leads to charge 
quantisation even with just a double barrier system, due 
to the insulating nature of the Luttinger wire in the pres- 
ence of any barrier. 

Although for the case of a rigidly sliding potential 
U{x — vt), where the periodicity in time T is related to 
the periodicity in space L by T = nL/v, n ^integer, 
it can be proved using topological arguments that the 
charge is quantised [1] , it is not so obvious that those ar- 
guments hold for more general cases where the Galilean 
principle does not hold. Several explicit examples were 
considered by Niu [2] but quantisation has been shown 
only for potentials which are spatially periodic. 

In this letter, we derive a simple formula for the 
pumped charge in terms of the reflection and transmis- 
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where r is the time period of the perturbation in the po- 
tential which causes the charge to be pumped and r, 6 
and t, (j) are the reflection and transmission amplitudes 
and phases respectively. We show that quantisation oc- 
curs whenever the contribution from the second term 
vanishes. 

We also explicitly demonstrate the quantisation of the 
pumped charge, due to interference effects, when we in- 
crease the number of tunnel barriers (n;,) or equivalcntly 
increase the number of dots (n^ — 71^ — 1) through which 
the current is measured. We work with open quantum 
dots with large transmissions and explicitly show how the 
charge pumped through the device changes with the num- 
ber of dots. It is expected [1,2] that for a fully spatially 
periodic system, (rif, or Ud — > oo), the charge pumped will 
be quantised. However, here we see that for reasonable 
values of the barrier strengths and pumping strengths, 
(almost) quantisation occurs even with 4-6 dots. We 
also see that the pumped charge, as a function of the 
number of dots, shows an interesting step-like behaviour. 
The pumped charge increases with the number of barri- 
ers and then saturates near an integer. Beyond that, it 
rises again with the number of barriers, till it saturates 
at the next integer. This tendency towards quantisation 
is also seen as a function of the amplitude of the pumping 
potential. This clearly indicates the special stability that 
occurs when an integer number of electrons are pumped 
through the system. Note that this stability is purely 
quantum mechanical in origin and is not due to interac- 
tions or CB physics, which allows electrons to be added 
only one by one to the dot. Here, we are completely ig- 
noring interaction effects and there is no CB since the 
dots are well-coupled to the leads. 

Although we use (5-function barriers for the explicit cal- 
culation, we expect the results to be robust to changing 
the form of the barriers. Following the work of Ref. [15], 
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wc also expect these results to be robust to weak disorder 
and to weak interactions. 
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FIG. 1. Schematic diagram of a multiple dot system (nj 
dots) defined on a two dimensional electron gas. The barriers 
forming the dots are denoted as Vi and the gate voltages 
controlling the density in the dots are denoted as Vg- 

We start with a system of coupled quantum dots as 
shown in Fig. 1. The barriers forming the dot are peri- 
odically modulated as 

Vi = Vi = Vo + Vp cos(cjt), i < ni,/2 for Ub = even, 

i < nf,/2 + I for rtb = odd, 
Vi =V2 = Vq + Vp cos{ojt + (5), for the remaining. (2) 

Here w is related to the time period as w = r/27r and 
S is the phase difference between the two time-varying 
potentials. Such a potential breaks the parity symmetry 
and allows the shape of all the dots to be varied. If 
the dots are coupled to the leads by a single channel 
quantum point contact, then it is sufficient to treat the 
dot within a one-dimensional effective Hamiltonian. The 
width of the dots (effectively the width of the quantum 
well that we use) is given by a. We are mainly interested 
in the region where Vq < Ep since we are in the resonant 
tunneling regime and not in the CB regime. 

The effective single channel S'-matrix for this system 
of rib barriers can be written as 



S = 



re 



(3) 



where the parameters r,t,r',9,9' and <j) are functions of 
the Fermi energy Ep and the amplitudes of the time- 
varying potentials Vi{t). Their explicit forms can be 
found, in terms of the parameters of a single well, (in 
the adiabatic limit), by solving the time- independent 
Schrodingcr equation for the potential Vi{t) given in Eq. 
2, for each value of t. The reflection amplitudes are not 
the same because the time-varying potentials explicitly 
violate parity. The potential also violates time-reversal 
invariance. But since in the adiabatic approximation, we 
are only interested in snapshots, at each value of the time, 
the Hamiltonian is time-reversal invariant and hence, the 
transmission amplitudes are the same for the 12 and 21 
elements in the S-matrix. 

By the Brouwer formula [6], the charge pumped can 
directly be computed from the parametric derivatives of 
the S-matrix. For a single channel, it is given by 
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where Sij denote the matrix elements of the S-matrix and 
V\ and V2 are the time derivatives of the V\ , V2 given in 
Eq. 2. For the form of the S-matrix given in Eq. 3, this 
is just 



Q = 



2lT 



(5) 



Thus, the pumped charge is directly related to the ampli- 
tudes and phases that appear in the scattering matrix. 
Note that Q can also be written in the form of Eq. 1 
where the first term is clearly quantised since e*^ has to 
return to itself at the end of the period. So the only pos- 
sible change in 6 can be in integral multiples of 2tt. The 
second term is the 'dissipative' term which prevents the 
perfect quantisation. It is easy to see the analogy of Eq. 
1 with Eq. 19 of Ref. [7]. 

The form in Eq. 5 also indicates that Q is quantised 
whenever either r or t is zero throughout the period. 
When the Fermi energy lies in a gap, {t = O.r = 1) 
the charge is quantised. This is what happens for spa- 
tially periodic potentials as discussed in Ref [1,2]. The 
charge is also quantised when there is total transmission 
(f = l,r = 0) through almost the whole period. This 
is essentially the case studied in Ref [12], where they 
find quantised charges whenever the pumping contour 
encloses almost all of the resonance line. 

In the rest of the letter, we compute the transmission 
and reflection coefficients, the phases and the quantised 
charges for various cases, and see that the pumped charge 
is almost quantised even when the number of barriers is 
quite small. We study the variation of Q function 
of the number of barriers, as a function of Ep and as a 
function of the pumping amplitudes. We also study Q as 
a function of the separation a between the barriers and 
as a function of the phase difference 5. 

Strictly speaking, to remain within adiabatic approxi- 
mation under which the Brouwer formula is derived, the 
energy level spacing in the dots A has to be larger than 
the energy scale defined by the frequency of the time- 
varying parameter Ei_j = Ulj. It is only under this ap- 
proximation that the snapshot picture of studying the 
static S'-matrix for different time points within the pe- 
riod is valid. A better approach to go beyond the adia- 
batic approximation [16,17] is to use the Floquet states. 
However, here we use the adiabatic approximation even 
in the continuum limit [Ep > Vq) where the energy levels 
are almost continuous and A ^ 0, and a few energy lev- 
els cross the Fermi level within a period. For sufSciently 
small u!, we expect this approximation to still yield qual- 
itatively correct results. 

• Single dot case or Ub = 2 : Here, we compute the scat- 
tering matrix for two ^-function barriers at a distance a 
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apart. Following Ref. [11], to obtain numerical results, we 
set a = 4 and w = 1. We find, however, that our results 
are independent of lo and hence w can be made as small as 
we wish. Our energy units are set hy Ti = 2m = fc^ = 1, 
where fcs is the Boltzmann constant. So for a — 100 A°, 
which is a typical value of the mean free path in GaAs, 
the energy unit in our system of units is E — b.GmeV, 
which corresponds to a temperature of T = 65°K. With 
these units, we set Vq — 1 and Vp — 0.4. We have also 
set the phase difference (5 = 7r/2 to obtain the maximum 
pumped charge. Here, we essentially reproduce the re- 
sults of Wei et al [11]. We have also checked that the 
peak in the small pumped current, occurs at the Fermi 
level, when one transmission maxima passes through the 
Fermi level in one period. Also, since the transmission 
does not become small (for a double-barrier system, the 
transmission does not fall to zero after reaching a max- 
imum), there is a large dissipative term in the pumped 
charge, which explains why the charge pumped is small. 
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FIG. 2. Q and the average transmission amplitude per pe- 
riod versus Ef for nt = 6 barriers. We have set Vo = 1, 
Vp = 0.4, = 1, a = 4 and 5 = tt/2. 

• Multiple dot case or ni, > 2 : Here, we have computed 
the S'-matrix, and obtained the transmission and reflec- 
tion coefficients, their phases and the pumped charge for 
Ub ranging from 3 to 14 ( 2 to 13 quantum dots). For each 
value of Tif,, the pumped charge and the average transmis- 
sion amplitude (obtained by integrating the transmission 
amplitude over one period and dividing by the period), 
is plotted as a function of the Fermi energy Ep, between 
Ep < Vq to E > Vq. This is shown in Fig. 2. for a typical 
case {rih — 6). The pumped charge is a maximum just 
when the transmission rises from zero. Since Ep <C Vq is 
the CB limit, which is not the limit we are studying, our 
main focus is on the charge pumped when E ^ Vq (the 
second peak in Fig. 2). Peaks at higher values of Ep 
become progressively smaller. However, we often work 
with the first peak also for illustrative purposes, since we 
have not included CB in our formalism. 

The magnitude of the (position of maximum as a func- 
tion of Ep) pumped charge as a function of rif, is plotted 
in Fig. 3 for the first and second peaks. Clearly, Q 
increases with the number of barriers. The first peak 
reaches quantisation with just — 6 whereas the sec- 
ond peak requires rib — 9. The most significant feature 




FIG. 3. Q (in units of e) as a function of the number of 
barriers nt- The top line is for the first peak and the line 
below is for the second peak. The parameter specifications 
are the same as for Fig. 2. 



here is the step-like structure or plateau structure near 
integer quantisation. (The quantisation for the first peak 
can be clearly seen, at 1 for Ub = 6, at 2 for Ub = 10 and 
at 3 for Ub = 14. For the second peak, we have only gone 
upto the value of one for Q, since it is progressively more 
cumbersome to go to larger number of barriers. But we 
have checked that the trend is the same.) For Ub oo, 
the system would be spatially periodic and perfect quan- 
tisation would have been expected. However, what is 
interesting and unexpected is the formation of plateaux 
for small values of Ub- 




FIG. 4. Transmission and reflection amplitudes {t and r) 
and phases (0 and 9) and the pumped charge Q as a function 
of time (for a single period r = 2tt) for nt = 10 barriers (left 
panel) and = 4 barriers (right panel). Q is measured in 
units of e snd (f) and 9 in radians. 

The correlation of the pumped charge with the values 
of the phases and amplitudes of the reflection and trans- 
mission coefRcients is shown in Fig. 4 for a couple of 
typical cases. In the top left panel, Q is plotted along 
with the 9 and (/> as a function of time for one full pe- 
riod, for the system with nt, = 10. For this system, we 
know from Fig. 3, that the total pumped charge is (al- 



3 



most) quantised at Q = 2. From the figure, we see that 
(j) comes back to itself at the end of one period, but 9 
has a discontinuity. It changes by 2 x 27r in one period. 
Thus, the first term in Eq. 1 gives a factor of 2. The 
criterion for quantisation is that the second term should 
vanish. The correlation of the quantisation with the (ap- 
proximate) vanishing of the second term can also be seen 
by looking at the figure. In the left panel, we note that 
where the transmission amplitude t is large, both 9 and 
(j) change very little and almost symmetrically. On the 
other hand, where 9 and (f) change very rapidly, the am- 
plitude t is very small. In the right panel, for contrast, 
we have studied nf, = 4 where Q is not quantised. We see 
that the change in 9 is 27r so that the first term in Eq. 
1 gives 1. But here, there are rapid and non-symmetric 
changes in 9 and (p when t is large. Hence, here, the sec- 
ond term is non-zero and there is no quantisation of Q 
near unity. Features similar to the left and right panels 
are consistently seen whenever there is quantisation and 
whenever there is no quantisation respectively. 

The qualitative features described above do not change 
when we change the ratio of to Vq- In fact, as Vp/Vo 
increases, we find that the value of Q and the tendency 
towards quantisation increases. This is seen in Fig. 5. 
We have also checked that Q is periodic in the separation 
a (shown as inset in Fig. 6.). In the weak pumping limit, 
- i.e., when the amplitude Vp is suffiiciently small, we 
expect the charge pumped to be proportional to sin S [6] , 
but as Vp increases, the sinusoidal shape is expected to be 
distorted. This feature is seen in Fig. 6. For Vp ~ 0.05, 
the dependence on S is sinusoidal, but here Q is quite 
small. As Vp increases, Q increases, but there is also an 
increasing distortion of the sinusoidal shape. 
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FIG. 5. Q as a function of the amplitude Vp of the 
time- varying potential for different ni,. 

To summarise, in this letter, we have shown that the 
pumped charge shows an interesting step-like behaviour 
as we increase the number of dots through which the 
charge is pumped and also as we increase the pumping 
amplitude. We have derived and demonstrated the rela- 
tion beween Q and the transmission and reflection coeffi- 
cients of the effective 5'-matrix. Q is quantised whenever 
the contribution from the transmission amplitude van- 
ishes - i.e., whenever there is no dissipation. The main 
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FIG. 6. Pumped charge as a function of the phase difler- 
ence 5 for n;, = 6 barriers. The inset shows it as a function of 
the separation a between barriers. 

point of this letter is to emphasize that this can occur due 
to quantum interference, even for transmission through 
a few (4-6) dots and not only for an infinite spatially 
periodic system. The experimentally testable prediction 
here is that the pumped charge increases as the number 
of barriers through which is is transmitted increases, and 
reaches quantisation with a few barriers. 
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